An artificial dynamic heterogeneity in a glass-former below Tg 
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We report the first large scale molecular dynamics simulations of the effect of the photoisomer- 
ization of probe molecules on the dynamics of a glassy molecular material. We show that the 
photoisomerization of the probe molecules induces the creation of an artificial dynamical hetero- 
geneity inside the matrix. We then study the modification of the dynamics and find, in agreement 
with recent experiments, an important increase of the diffusion leading to liquid like diffusion below 
the glass transition temperature. We find that the photoisomerization process controls the hetero- 
geneity and the non-Gaussian parameter of the material, leading to extremely rapid variations of 
these quantities. 
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A central question in the field of the glass-transition 
is the relation between dynamic heterogeneities and the 
relaxation dynamics of supercooled liquids and glassy 
materials [TJ [51 [3]. In order to shed some light on 
this relation, we investigate various methods to modify 
the dynamic heterogeneities, and then to observe the 
resulting effects in the relaxation dynamics. Because it 
induces a cutoff in the correlation length of any coopera- 
tive mechanism present in the material, confinement has 
been seen as a promising direction of study. The liquid 
may be confined inside pores, micelles, or nanoparticles 
may be dispersed inside the liquid. However in most 
cases the surface effects predominate in the origin of the 
dynamical properties of confined supercooled liquids. 
For example, depending on the surface properties of 
the confining medium the dynamics is slowed down or 
accelerated [4j ^ ,7, . 

In this letter we present a different and more di- 
rect method to modify the dynamic heterogeneities 
inside a supercooled or a glassy material. We show that 
the photoisomerization of probe molecules inside the 
glassy or supercooled matrix induces the creation of 
an artificial dynamical heterogeneity inside the matrix. 
We then study the modification of the dynamics and 
find, in accordance with recent experimental results[8j, 
an important increase of the diffusion. We find a 
strong modification of the plateau of the mean square 
displacement (MSD), suggesting induced cage breaking 
process as a cause for the acceleration of the dynamics. 
We find that the photoisomerization process controls the 
heterogeneity and the Non-Gaussian parameter (NGP) 
of the material, leading to extremely rapid variations of 
these quantities. We simulate the photoisomerization of 
dispersed red {DRl) molecules (probes) inside a matrix 



of methylmethacrylate {MM A) molecules {Cc,02Hg). 
We study the effect of the probe isomerization on the 
matrix dynamics and on the dynamical heterogeneity 
inside the matrix. Under continuous holographic expo- 
sure of thin polymer films doped with azo-dyes Rochon 
et al. and Kim et al.[ni HO] observed the formation of 
surface relief gratings (SRG) corresponding to massive 
photoinduced mass transport [HI [TTJ [121 ESI well be- 
low the glass transition temperature of the material. 
Various theories have been proposed to explain this 
effect [H [m [13 [m [ig. Recently SRG have also been 
observed in molecular glasses [13 HOI [H] suggesting 
universality in the mechanism of SRG formation. An 
increase of the local diffusion has been suggested [lH [17] 
as a possible cause for SRG formation but is still a 
matter of controversy. 

The present simulations were carried out for a sys- 
tem of 700 MM A molecules and 1 or 7 DRl molecules. 
Our simulations use the Gear algorithm with the quater- 
nion method [Hj to solve the equations of motions with 
the following potential functions [121 IMj- The time step 
was chosen equal to 0.5 10~^^s. The density was set con- 
stant in our simulations at \.\%g / cm? . The temperature 
was controlled using a Berendsen thermostat ^5]. The 
observed maximum fluctuations were of +/— SK. The 
DRl molecules are in some simulations treated as rigid 
bodies in the trans or cis isomer in order to study the 
dynamics without photoisomerization. In other simula- 
tions, in order to study the effect of photoisomerization, 
the DRl molecules change periodically and continuously 
from the trans to the cis isomer and then from the cis to 
the trans isomer. In this case the DRl molecule changes 
from trans to cis uniformly in a time ti = 300/s. This 
value is typical of the isomerization process |26j . Then it 
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stays in the cis configuration during a time ^2, changes 
uniformly from cis to trans in a time ti and stays in the 
trans configuration during a time ^2- Results presented 
in this letter correspond to t2 = ti then a short waiting 
time value. We have then tested the effect of the time 
variation of ti and t2 on the matrix dynamics |27j . 
We find [27^ that the dependence is roughly linear with 
time variations. Typical experimental light exposure 
leads to a value of 1ms between isomerizations. Taking 
also into account the weak DRl concentration of our 
simulations, we estimate that our diffusion coefficient 
must be multiplied by a rough factor 10^^ if compared 
with typical SRG experiments. We define the mobility 
Mj,to(0 of the atom i at time to with a characteristic 
time t, as the displacement of the atom i during the 
time ^pg]: 

H.^toit) = \r^{t + ta) - r,{ta)\ [1] 

The dynamical heterogeneity are characterized [29l [30j 
[3TJ |32] by the transient aggregation of the molecules of 
highest mobility in the liquid. This concerns typically the 
5 percent most mobile molecules. Similarly a transient 
aggregation of the molecules of low mobility is observed 
[28l[33] but with different characteristic times. These ag- 
gregations are the fingerprints of the cooperativity of mo- 
tions inside the supercooled hquid and disappear around 
the melting temperature [221 133 1311 132 • In order to 
quantify the aggregation of the most mobile molecules we 
define the function: 

rRa 

I^it)= / (.g™-™,t(r)/.g(r) - l)Wdr [2] 
Jo 

and similarly for the least mobile molecules: 

I~{t)= / {9im^im.t{r)/9{r)-l)'^TTr^dr [3] 

Where gm~m,t{f) is the radial distribution function 
(RDF) between the most mobile matrix molecules. 
And gim-im,t{'r') is the RDF between the least mobile 
molecules, while g{r) is the RDF between random matrix 
molecules. The time t appearing in these functions is the 
time used in the definition (1) of the mobility iii{t). This 
time is then used in the selection of the most or the least 
mobile molecules. The cut-off Rc = 6 A has been chosen 
to reduce noise effect. For text clarity we will name 
functions I^^^(t) intensity of the aggregation. When 
the time t corresponds to the phenomena of physical 
importance for the motion of molecules, the aggregation 
appears. The aggregation of mobile molecules then 
appears usually [29] for t = t* which corresponds to 
the maximum value of the NGP. This time scale corre- 
sponds roughly to the mean time used by a molecule to 
escape the cage created by its neighbors. For the least 
mobile molecules the characteristic time is usually larger 
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FIG. 1: Intensity of the mobiles heterogeneities with = 6 
A, and of the slow heterogeneities, for the centers of masses of 
the matrix molecules. The photoisomerization is switched on. 
The temperature is 180 K, then well below the glass transition 
temperature (Tg — 400 A"). The simulation box contains one 
DRl molecule and 700 MM A molecules. 

and this difference increases when the temperature 
decreases. This effect has been related to the breaking 
of the Stokes Einstein relation [TJ |34]. These behavior 
have been observed in a number of supercooled liquids 
[H m [301 [311 [32]. When the photo-isomerization is 
switched off, we observe this behavior in the MM A + 
DRl material investigated here [27] . 

When the photo-isomerization is switched on, the 
behavior of the dynamic heterogeneities changes totally. 
Figure 1 shows the intensity versus time I~^{t) (eq. 2) of 
the aggregation of the most mobile matrix molecules in 
this case . The temperature of this simulation (180 AT) is 
largely below the glass-transition temperature. At this 
temperature the time scale of the heterogeneity when the 
isomerization is switched off is larger than Ins [27]. The 
time scales in figure 1 are much shorter. The intensity 
I'^it) oscillates rapidly between maxima and minima. 
The most mobile matrix molecules aggregate when the 
DRl probe molecule undergoes a trans-cis isomeriza- 
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tion. And these aggregates disappear partly when the 
probe molecule undergoes a cis-trans isomerization. To 
our knowledge the time scale of the heterogeneity ob- 
served here is much shorter than any time scale reported 
to date. This time scale corresponds exactly to the time 
scale introduced in the simulation for the isomerization. 
This result shows that the heterogeneity in the matrix 
is here controlled by the isomerization of the probe 
molecule. Figure 1 also shows the intensity I^{t) of 
the aggregation of the least mobile molecules. I~{t) 
displays peaks at periodic times. The slow heterogeneity 
is maximal when the DRl is in the trans isomer. The 
slow heterogeneity then disappears completely when the 
DRl undergoes a trans-cis isomerization. We observe in 
figure 1 that the intensities /+ {t) and (t) , correspond- 
ing to the aggregation of mobile and slow molecules, 
follow an opposite behavior. When the DRl undergoes 
a cis-trans isomerization the slow heterogeneity appears 
and increases sharply to a maximum value. At the same 
time the mobile heterogeneity decreases sharply. Then 
when the DRl undergoes a trans-cis isomerization 
the opposite behavior appears. The slow heterogeneity 
collapses and the mobile heterogeneity increases to a 
maximum value. We observe here for the first time, 
a total time separation between the mobile and slow 
heterogeneity. The time evolution of the heterogeneity 
coincides with the time evolution of the isomerization 
process. When the photo-isomcrization is switched off 
the heterogeneous behavior of this material at this tem- 
perature corresponds to characteristic times larger than 
Ins. No heterogeneous behavior is then observed on 
the time scales of figure 1 when the photo-isomerization 
process is switched off. In the material studied here, the 
photo-isomerization of the DRl then controls totally 
the slow and mobile heterogeneities. These results show 
the possibility of controlling the heterogeneity inside a 
material with an external input. 



We observe in Figure 2, peaks and dips superimposed 
on a main smooth curve. The main smooth curve has 
an usual shape for a2{t). But the time scale correspond 
to a viscous liquid at a much higher temperature. The 
characteristic time that corresponds to the maximum of 
the main smooth curve is t\ — 8 ps. While at the same 
temperature the characteristic time t* of a2{t) is higher 
than Ins when the isomerization is switched off. The 
characteristic time of the main curve ti is then very un- 
usual. Superimposed on this main curve, peaks and dips 
appear with the characteristic times of the trans and cis 
lifetimes used in our simulation. The time scales of these 
peaks and dips are very short. The correspondence be- 
tween these time scales and the trans and cis life times 
show that this behavior is also controlled by the isomer- 
ization process. The DRl isomerization controls then 
the NGP value, which is an important characteristic of 
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FIG. 2: Non Gaussian parameter a2(t) of the centers of 
masses of the matrix molecules when the photoisomerization 
is switched on. The temperature is 180K. The simulation box 
contains one DRl molecule and 700 MM A molecules. 



the non-Markov dynamics. Figure 3 shows the MSD of 
the matrix molecules when the DRl undergoes or not 
photo-isomerization. At this temperature when there is 
no isomerization, the motion is stopped on the time scale 
of the simulation. This is shown on the lower curve of 
Figure 3. The plateau of the MSD begins around Ips and 
continues up to the larger times of the simulation. How- 
ever when the isomerization is switched on, the plateau 
of the MSD disappears and is replaced by oscillations. 
Figure 3 shows that for large times the MSD follows the 
Einstein law |35j for Brownian motion (dashed line) . For 
large times we then observe a liquid- like diffusion law. 
The diffusion coefficients D calculated from this law are 
displayed in Figure 4. These diffusion coefficients are 
typical of viscous liquids and not of glassy materials. We 
then observe in Figure 3 a liquid-like diffusion mechanism 
at a temperature below Tg. This strange comportment 
is validated by recent experimental data [8]. The oscil- 
lations in the MSD correspond to the cis-trans isomer- 
ization process. The MSD is maximum for the trans- 
cis isomerization and minimum for the cis-trans isomer- 
ization. The comportment of the matrix for trans-cis 
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FIG. 3: Mean square displacement < r'^{t) > versus time 
at a temperature of 180K, when the photoisomerization is 
switched on (upper curve) or ofT (lower curve) . The simulation 
box contains 7 DRl molecule and 700 MM A molecules. 



isomerization approaches the comportment of a liquid. 
The NGP decreases showing that the dynamics tends to 
a Markov dynamics; The slow heterogeneity disappears 
and the diffusion coefficient D approaches viscous liquid 
values. However the presence of mobile heterogeneity 
shows that the matrix does not behave totally as a nor- 
mal liquid. Amazingly, the comportment of the matrix 
for cis-trans isomerization is also not the comportment 
of a glassy material. Approaching the comportment of 
a glass, the slow heterogeneity increases, the NGP in- 
creases and the MSD decreases. However the MSD is 
higher than the plateau of the glass, as seen in figure 
3. Another striking difference is that the characteristic 
time of the slow heterogeneity is not increasing with tem- 
perature decrease. The characteristic times for slow and 
mobile heterogeneity and NGP are here imposed by the 
isomerization characteristic times. When the isomeriza- 
tion trans-cis is switched on, it imposes a cage breaking 
process. Figure 3 shows that the plateau, which is a char- 
acteristic of the cage effect, disappears. This cage break- 
ing process follows then the characteristic times of the 
isomerization. When the isomerization is switched off. 
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FIG. 4: Diffusion coefficient for the matrix, as a func- 
tion of the inverse of temperature, when the photoisomer- 
ization is switched on: Empty circles; or off: Full circles. 
The dashed line fit (photoisomerization switched off) corre- 
sponds to an Arrhenius dependence D — D\exp[—Ea\ /T) 
with Eai ~ iOOOK. When the isomerization mechanism 
is switched on, the diffusion coefficient increases. The dot- 
ted line fit (photoisomerization switched on) follows a law 
of the form: D = Diexp{-Eai/T) + D2exp{~Ea2/T) with 
Ea2 = 250K. This result suggests that the isomerization 
adds a new Arrhenius process to the original dynamics of the 
matrix. We estimate that D must be multiplied by a rough 
factor fO~^ if compared with typical SRG experimental con- 
ditions. 

the diffusion coefficient of the matrix molecules, doped 
with trans or cis DRl isomers, are very similar j27j. The 
trans DRl isomer has been chosen in the figures. 

Figure 4 shows the diffusion coefficient of the matrix 
molecules versus temperature when the isomerization is 
switched on or off. When the isomerization is switched 
off we have an Arrhenius evolution of the diffusion 
coefficient. We observe a typical Arrhenius dependence 
D = Diexp{~Eai/T) with Eai = 4000K. The diffusion 
coefficients in Figure 4 may also be fitted with an 
MCT law: D = Do{T ~ Tc)"' with Tc = 447.2^ 
and 7 = 2.3. The MCT fit leads to an extrapola- 
tion of Tg around iOOK. When the isomerization 
mechanism is switched on, the diffusion coefficient 
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increases. We obtain the sum of two Arrhenius mech- 
anisms. The diffusion coefficient follows a law of the 
form: D = Diexp{-Eai/T) + D2exp{-Ea2/T) with 
Ea2 = 2bQK. This result suggests that the isomerization 
adds a new Arrhenius process to the original dynamics 
of the matrix. 

In conclusion, our study shows that the isomeriza- 
tion of probe molecules in glassy materials creates 
a dynamical heterogeneity inside the material. This 
heterogeneity is controlled by the isomerization process. 
This result suggests a new method to study the effect 
of the heterogeneity on the material dynamics, using 
the photoisomerization of probe molecules to control 
the heterogeneity. To our knowledge, this is a unique 
case where the slow heterogeneity and the mobile 
heterogeneity appear at different times and thus are not 
present together. This result shows that the two kinds 
of heterogeneity are not directly related, i.e. that the 
slow heterogeneity is not the cause of the presence of 
the mobile heterogeneity in supercooled liquids. These 
results suggest a possible control of the dynamical prop- 
erties of amorphous materials and of the glass-transition 
mechanism via the creation of heterogeneity which may 
lead to a new route for cryopreservation of biomaterials. 
The huge effect observed on the dynamical properties 
of the material suggests a strong relation between the 
presence of dynamical heterogeneity and the dynamical 
properties. This study also shows that both kinds of 
aggregations (slow and mobile) have a different effect 
on the dynamics. This result has to be taken into 
account in the Adam et GibbS jSj interpretations of the 
glass-transition when mobile heterogeneity are used as a 
measure of the configurational entropy. 
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